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Abstract 

This paper investigates packing and covering properties of codes with the rank metric. First, we 
investigate packing properties of rank metric codes. Then, we study sphere covering properties of rank 
metric codes, derive bounds on their parameters, and investigate their asymptotic covering properties. 
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Although the rank has long been known to be a metric implicitly and explicitly (see, for example, [1]), 
the rank metric was first considered for error control codes by Delsarte [2]. The potential applications of 
rank metric codes to wireless communications [3], public-key cryptosystems [4], and storage equipments 
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[5], [6] have motivated a steady stream of works [2], [5]— [19], described below, that focus on their 
properties. 

XI 

The majority [2], [5]-[7], [12], [13], [15], [17], [18] of previous works focus on rank distance properties, 
code construction, and efficient decoding of rank metric codes. Some previous works focus on the packing 
and covering properties of rank mettic codes. Both packing and covering properties are significant for 
error control codes, and packing and covering radii are basic geometric parameters of a code, important 
in several respects [20]. For instance, the covering radius can be viewed as a measure of performance: 
if the code is used for error correction, then the covering radius is the maximum weight of a correctable 
error vector [21]; if the code is used for data compression, then the covering radius is a measure of the 
maximum distortion [21]. The Hamming packing and covering radii of error control codes have been 
extensively studied (see, for example, [22], [23]), whereas the rank packing and covering radii have 
received relatively little attention. It was shown that nontrivial perfect rank metric codes do not exist in 
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[8], [9], [16]. In [10], a sphere covering bound for rank metric codes was introduced. Generalizing the 
concept of rank covering radius, the multi-covering radii of codes with the rank metric were defined in 
[11]. Bounds on the volume of balls with rank radii were also derived [19]. 

In this paper, we investigate packing and covering properties of rank metric codes. The main contri- 
butions of this paper are: 

• In Section [Till we study the packing properties of rank metric codes. 

• In Section llVl we establish further properties of elementary linear subspaces [18], and investigate 
properties of balls with rank radii. In particular, we derive both upper and lower bounds on the 
volume of balls with given rank radii, and our bounds are tighter than their respective counterparts 
in [19]. 

• In Section [V] we first derive both upper and lower bounds on the minimal cardinality of a code 
with given length and rank covering radius. Our new bounds are tighter than the bounds introduced 
in [10]. We also establish additional sphere covering properties for linear rank metric codes, and 
prove that some classes of rank metric codes have maximal covering radius. Finally, we establish 
the asymptotic minimum code rate for a code with given relative covering radius. 

II. Preliminaries 

A. Rank metric 

Consider an n-dimensional vector x = {xq, x%, . . . , x n _i) G GF(q m ) n . The field GF(q m ) may be 
viewed as an m-dimensional vector space over GF(g). The rank weight of x, denoted as rk(x), is 
defined to be the maximum number of coordinates in x that are linearly independent over GF(g) [7]. 
Note that all ranks are with respect to GF(g) unless otherwise specified in this paper. The coordinates 
of x thus span a linear subspace of GF{q m ), denoted as <5(x) or &(xq,xi, . . . ,x n -\), with dimension 
equal to rk(x). For any basis B m of GF{q m ) over GF(g), each coordinate of x can be expanded to an 
m-dimensional column vector over GF(g) with respect to B m . The rank weight of x is hence the rank of 
the mxn matrix over GF(g) obtained by expanding all the coordinates of x. For all x, y G GF(q m ) n , it 

def 

is easily verified that d R (x, y) = rk(x — y) is a metric over GF(q m ) n [7], referred to as the rank metric 
henceforth. The minimum rank distance of a code C, denoted as d R (C), is simply the minimum rank 
distance over all possible pairs of distinct codewords. When there is no ambiguity about C, we denote 
the minimum rank distance as d R . 

Both the matrix form [2], [6] and the vector form [7] for rank metric codes have been considered in 
the literature. Following [7], in this paper the vector form over GF(q m ) is used for rank metric codes 
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although their rank weight can be defined by their corresponding m x n code matrices over GF(g) [7]. 
The vector form is chosen in this paper since our results and their derivations for rank metric codes can 
be readily related to their counterparts for Hamming metric codes. 

B. Sphere packing and sphere covering 

The sphere packing problem we consider is as follows: given a finite field GF{q m ), length n, and 
radius r, what is the maximum number of non-intersecting balls with radius r that can be packed into 
GF(q m ) n ? The sphere packing problem is equivalent to finding the maximum cardinality of a code 
over GF(g m ) with length n and minimum distance d > 2r + 1: the spheres of radius r centered at the 
codewords of such a code do not intersect one another. Furthermore, when these non-intersecting spheres 
centered at all codewords cover the whole space, the code is called a perfect code. 

The covering radius p of a code C with length n over GF{q m ) is defined to be the smallest integer p 
such that all vectors in the space GF(q m ) n are within distance p of some codeword of C [23]. It is the 
maximal distance from any vector in GF{q rn ) n to the code C. That is, p = max xg (j F ( g ™), 1 {d(x, C)}. 
Also, if C C C , then the covering radius of C is no less than the minimum distance of C. Finally, 
a code C with length n and minimum distance d is called a maximal code if there does not exist any 
code C with the same length and minimum distance such that C C C. A maximal code has covering 
radius p < d — 1. The sphere covering problem for the rank metric can be stated as follows: given an 
extension field GF{q m ), length n, and radius p, we want to determine the minimum number of balls of 
rank radius p which cover GF(q m ) n entirely. The sphere covering problem is equivalent to finding the 
minimum cardinality of a code over GF(q m ) with length n and rank covering radius p. 

III. Packing properties of rank metric codes 

It can be shown that the cardinality K of a code C over GF(g m ) with length n and minimum 
rank distance d R satisfies K < min jg m (™- rf R+ 1 ) ; q,n(m-d R +i) j_ re f er t bound as the Singleton 
bounc^ for codes with the rank metric, and refer to codes that attain the Singleton bound as maximum 
rank distance (MRD) codes. 

For any given parameter set n, m, and d R , explicit construction for linear or nonlinear MRD codes 
exists. For n < m and d R < n, generalized Gabidulin codes [13] can be constructed. For n > m and 
d R < m, an MRD code can be constructed by transposing a generalized Gabidulin code of length m 

'The Singleton bound in [6] has a different form since array codes are defined over base fields. 
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and minimum rank distance d R over GF(q n ), although this code is not necessarily linear over GF(g m ). 
When n = Im (I > 2), linear MRD codes of length n and minimum distance d R can be constructed by a 
cartesian product Q l of an (m, k) generalized Gabidulin code Q. Although maximum distance separable 
codes, which attain the Singleton bound for the Hamming metric, exist only for limited block length over 
any given field, MRD codes can be constructed for any block length n and minimum rank distance d R 
over arbitrary fields GF(q m ). This has significant impact on the packing properties of rank metric codes 
as explained below. 

For the Hamming metric, although nontrivial perfect codes do exist, the optimal solution to the sphere 
packing problem is not known for all the parameter sets [22]. In contrast, for rank metric codes, although 
nontrivial perfect rank metric codes do not exist [8], [9], MRD codes provide an optimal solution to the 
sphere packing problem for any set of parameters. For given n, m, and r, let us denote the maximum 
cardinality among rank metric codes over GF{q m ) with length n and minimum distance d R = 2r + 1 as 
A R (q m ,n,d R ). Thus, for d R > min{n, m}, A R (q m ,n,d R ) = 1 and for d R < min{n, m}, A R (q m ,n,d R ) = 
mm jgm.(n-d R +i) ^n(m-d R +\) j j^ote t ^ a i the maximal cardinality is achieved by MRD codes for all 
parameter sets. Hence, MRD codes admit the optimal solutions to the sphere packing problem for rank 
metric codes. 

IV. Technical results 
A. Further properties of elementary linear subspaces 

The concept of elementary linear subspace was introduced in our previous work [18]. It has similar 
properties to those of a set of coordinates, and as such has served as a useful tool in our derivation of 
properties of Gabidulin codes (see [18]). Although our results may be derived without the concept of 
ELS, we have adopted it in this paper since it enables readers to easily relate our approach and results 
to their counterparts for Hamming metric codes. 

If there exists a basis set B of vectors in GF(q) n for a linear subspace V C GY{q m ) n , we say V is 
an elementary linear subspace and B is an elementary basis of V. We denote the set of all ELS's of 
GF(g m ) n with dimension v as E v (q m ,n). The properties of an ELS are summarized as follows [18]. A 
vector has rank < r if and only if it belongs to some ELS with dimension r. For any V G E v (q m ,n), 
there exists V £ E n _ v (q m ,n) such that V © V = GF(q m ) n , where © denotes the direct sum of two 
subspaces. For any vector x G GF(q m ) n , we denote the projection of x on V along V as xy, and we 
remark that x = xy + x v . 
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In order to simplify notations, we shall occasionally denote the vector space GF(q m ) n as F. We denote 
the number of vectors of rank u (0 < u < min{m,n}) in GF(q m ) n as N u (q m ,n). It can be shown that 
N u (q m ,n) = [™]a(m,u) [7], where a(m,0) d = 1 and a(m,u) = Ui=o(l m ~ ?*) for u > L The ["] 
term is often referred to as a Gaussian polynomial [24], defined as [™] = a(n,u)/a(u,u). Note that 
\E u (q m ,n) \ = [™] does not depend on m. 

Lemma 1: Any vector x G GF(q m ) n with rank r belongs to a unique ELS V G E r (q m ,n). 

Proof: The existence of V € E r (q m , n) has been proved in [18]. Thus we only prove the uniqueness 
of V, with elementary basis {vj}[~Q, where Vj G GF(g) n for all i. Suppose x also belongs to W, 
where W G E r (q m ,n) has an elementary basis {wj}^~Q, where Wj E GF(q) n for all j. Therefore, 
x = Yll=o a i v i = Sj=o^' w i' where aj,foj G GF(g m ) for < i,j < r — 1. By definition, we have 
S(x) = S(ao, . . . , a r _i) = S(&o, ■ • • , K-i), therefore bj's can be expressed as linear combinations of 
en's, i.e., bj = J2iZo c j,i a i where c jti G GF(q). Hence x = ^1=0°^' where u« = Y^Zq c j,i w j for 
< i < r — 1 form an elementary basis of W. Considering the matrix obtained by expanding the 
coordinates of x with respect to the basis {ai}™^ 1 , we obtain Vj = Uj, and hence V = W. ■ 

Lemma Q] shows that an ELS is analogous to a subset of coordinates since a vector x with Hamming 
weight r belongs to a unique subset of r coordinates, often referred to as the support of x. 

In [18], it was shown that an ELS always has a complementary elementary linear subspace. The 
following lemma enumerates such complementary ELS's. 

Lemma 2: Suppose V G E v (q m ,n) and A C V is an ELS with dimension a, then there are q a ( v ~ a ^ 
ELS's B such that A B = V. Furthermore, there are q a \"- a > [ v ] such ordered pairs (A, B). 

Proof: First, remark that dim(23) = v — a. The total number of sets of v — a linearly independent 
vectors over GF(g) in V\A is given by N = (q v -q a )(q v -q a+1 ) ■ ■ ■ {q v -q v ~ l ) = q a{v ~ a) a{v-a,v-a). 
Note that each set of linearly independent vectors over GF(g) constitutes an elementary basis set. Thus, 
the number of possible B is given by N divided by a(v — a,v — a), the number of elementary basis sets for 
each B. Therefore, once A is fixed, there are q a ( v ~ a > choices for B. Since the number of a-dimensional 
subspaces A in V is ["] , the total number of ordered pairs (A, B) is hence q a ( v ~ a ) . ■ 

Puncturing a vector with full Hamming weight results in another vector with full Hamming weight. 
Lemma [3] below shows that the situation for vectors with full rank is similar. 

Lemma 3: Suppose V G E v (q m ,n) and u G V has rank v, then rk(u_4) = a and rk(us) = v — a for 
any A G E a {q m , n) and B G E v _ a (q m , n) such that A@B = V. 

Proof: First, 114 G A and hence rk(u^) < a by [18, Proposition 2]; similarly, rk(ug) < v — a. 
Now suppose rk(u^) < a or rk(ug) < v — a, then v = rk(u) < rk(u_4) + rk(ug) < a + v — a = v. ■ 
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It was shown in [18] that the projection 114 of a vector u on an ELS A depends on both A and 
its complement B. The following lemma further clarifies the relation: changing B always modifies U4, 
provided that u has full rank. 

Lemma 4: Suppose V G E v (q m ,n) and u € V has rank v. For any A G E a (q m ,n) and B G 
E v _ a (q m ,n) such that A B = V, define the functions f u (A,B) = 114 and g u (A,B) = u#. Then 
both / u and g u are injective. 

Proof: Consider another pair (A',B') with dimensions o and v — a respectively. Suppose A' 7^ A, 
then u_4' 7^ 114. Otherwise 114 belongs to two distinct ELS's with dimension a, which contradicts 
Lemma Q] Hence 114/ 7^ u_4 and uq> = u — 114- 7^ u — 114 = u#. The argument is similar if B' 7^ B. ■ 

B. Properties of balls with rank radii 

We refer to all vectors in GF(q m ) n within rank distance r of x € GF(q m ) n as a ball of rank radius r 
centered at x, and denote it as -B r (x). Its volume, which does not depend on x, is denoted as V r (q m ,n) = 
Ylu=o N u (q m , n). When there is no ambiguity about the vector space, we denote V r (q m ,n) as v(r). 

Lemma 5: For < r < min{n,m}, q r ( m + n ~ r ) < V r (q m ,n) < K^ 1 q r ^ m+n ~ r \ where K q A = 

nr=i(i-<r J ) [i8]. 

Proof: The upper bound was derived in [18, Lemma 13], and it suffices to prove the lower bound. 
Without loss of generality, we assume that the center of the ball is 0. We now prove the lower bound 
by constructing q r ( m + n - r ) vectors z G GF(q m ) n of rank at most r. Let x G GF(q m ) r let a subspace X 
of GF(q m ) such that dim(T) = r and G(x) C T. We consider the vectors y G GF(q m ) n ~ r such that 
S(y) C T. There are q mr choices for x and, for a given x, g r '( n ~ r ) choices for y. Thus the total number 
of vectors z = (x, y) G GF(q m ) n is q r ( m + n - r ) f an d since 6(z) C T, we have rk(z) < r. ■ 
We remark that both bounds in Lemma [5] are tighter than their respective counterparts in [19, Proposi- 
tion 1]. More importantly, the two bounds in Lemma [5] differ only by a factor of K q , and thus they not 
only provide a good approximation of V r (q m ,n), but also accurately describe the asymptotic behavior 
of V r {q m ,n). 

The diameter of a set is defined to be the maximum distance between any pair of elements in the set 
[22, p. 172]. For a binary vector space GF(2) n and a given diameter 2r < n, Kleitman [25] proved that 
balls with Hamming radius r maximize the cardinality of a set with a given diameter. However, when 
the underlying field for the vector space is not GF(2), the result is not necessarily valid [23, p. 40]. We 
show below that balls with rank radii do not maximize the cardinality of a set with a given diameter. 
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Proposition 1: For 2 < 2r < n < to, any 5 E E2 r {q m ,n) has diameter 2r and cardinality |5| > 

V r (q m ,n). 

Proof: Any S E E2 T {q m , n) has diameter 2r and cardinality q 2mr . For r = 1, we have Vi(g m , n) = 
1 + ^""(g-i""^ < 9 2m For ^ > 2, we have K-(g m , n) < K~ l q r(n+™)-r 2 by L emma|5] Since r 2 > 2 > 
- log g Kg, we obtain V r (q m , n) < «f < \S\. ■ 

The intersection of balls with Hamming radii has been studied in [23, Chapter 2], and below we 
investigate the intersection of balls with rank radii. 

Lemma 6: If < r, s < n and ci,C2 E GF(q m ) n , then |B r (cj) n5 s (c2)| depends on ci and C2 only 
through d R (ci,c 2 ). 

Proof: This follows from the fact that matrices in GF(g) mxn together with the rank metric form 
an association scheme [2], [26]. ■ 
Proposition 2: If < r, s < n, ci^jc'^c^ E GF(q m ) n and d R (ci,c 2 ) > ^(c^cQ, then \B r (ci) n 

a 5 (c 2 )| < |s P (ci)ns a (c^)|. 

Proof: It suffices to prove the claim when d R (cx,C2) = ^(c^c^) + 1 = e + 1. By Lemma |6j 
we can assume without loss of generality that ci = = 0, c 2 = (0, ci, . . . , c e , 0, . . . , 0) and c 2 = 
(co, ci, . . . , c e , 0, . . . , 0), where co, • • • , c e E GF(g m ) are linearly independent. 

We will show that an injective mapping <fi from B r (ci)r\B s (c2) to Br^'-^^B s {c' 2 ) can be constructed. 
We consider vectors z = (zo, Z\, . . . , Zn-i) E B r (ci) n B s (c 2 ). We thus have rk(z) < r and rk(u) < s, 
where u = (it , iti, . . . ,« n _i) = z-c 2 = (z — cq, — ci, . . . , z„-i). We also define z = (zi, . . . , z n -i) 
and u = (tti, . . . , u n -i). We consider three cases for the mapping <j>, depending on z and u. 

def 

• Case I: rk(u) < s — 1. In this case, $>(z) = z. We remark that rk(z — c 2 ) < rk(u) + 1 < s and 
hence 0(z) E B r (c9 n B s (c 2 ). 

def 

• Case II: rk(u) = s and rk(z) < r — 1. In this case, <p(z) = (zq — Co, Zi, • • • , z n -i). We have 
rk(0(z)) < rk(z) + 1 < r and rk ((/>(z) — c 2 ) = rk (z — c 2 ) < s, and hence 4>(z) E B r (cJ)nB s (c' 2 ). 

• Case III: rk(u) = s and rk(z) = r. Since rk(u) = s, we have z$ — cq E S(u). Similarly, since 
rk(z) = r, we have zo E <5(z). Denote dim((5(u, z)) as d (ci > s). For d > s, let ao, • • • , ad-i be 
a basis of <5(u, z) such that ao, . . . , a s -i E <5(u) and a s , ■ ■ ■ , ad-i £ This basis is fixed for 
all vectors z having the same z, i.e., it is fixed for all values of Zq. Note that co E <5(u, z), and may 
therefore be uniquely expressed as co = c u + c 2 , where c u E S(ao, ■ ■ ■ , cn s -i) = <5(u) and c 2 E 

def 

<5(a s , . . . , a^-i) C 6(z). If d = s, then c z = E 6(z). In this case, 0(z) = (zq — c z , zi,..., z n -i)- 
Remark that zq — c z E S(z) and hence rk(0(z)) = r. Also, z$ — c z = zq — cq + c u E <5(u) and 
hence rk(^>(z) — c' 2 ) = s. Therefore 0(z) E ^(c^) n 5 s (c 2 ). 
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It can be easily verified that <f> is injective, hence |J3 r (ci) n B s {c2)\ < | S r (cj_ ) n B s (c 2 )|. ■ 

Corollary 1: If < r, s < n, ci,C2,Cj,C2 G GF(g m )™ and d R (ci,c 2 ) > d R (c^,c 2 ), then \B r {c\) U 
B s {c 2 )\ > \B r (c[)UB s (c' 2 )\. 

Proof: The result follows from ^(ci) UB s (c 2 )| = v(r) + v(s) - \B r (ci) n B s (c 2 )\. ■ 

We now quantify the volume of the intersection of two balls with rank radii for some special cases, 
which will be used in Section IV-BI 

Proposition 3: If ci,c 2 G GF(q m ) n and d R (ci, c 2 ) = r, then |S r (ci) n #i(c 2 )| = 1 + (g m - q r ) [[] + 
(<f-l)[i]- 

Proof: By Lemma [T] the vector ci belongs to a unique ELS V G E r (q m ,n). First of all, it is easy 
to check that y = G -B r (ci) n -E>i(0). We now consider a nonzero vector y G -Bi(O) with rank 1. We 
have d R (y, ci) = r + 1 if and only if y ^ V and <5(y) ^ ©(ci). There are ^ q ~ 9 g ^ ~ g - such vectors. 
Thus, |B r ( Cl ) n Bi(c 2 )\ = 1 + N!(<r,n) - togf^) = 1 + (<f» - g r) [r] + {q r _ 1} ^ ■ 
Proposition 4: If ci,c 2 G GF(g m ) n and d R (ci,c 2 ) = r, then \B s (c x ) nB r _ 8 (c 2 )| = g s(r_s) Q for 
< s < r. 

Proof: By Lemma [6l we can assume that ci = 0, and hence rk(c 2 ) = r. By Lemma [Q c 2 belongs 
to a unique ELS V G E r (q m ,n). We first prove that all vectors y G -B s (0) n B r _ s (c 2 ) are in V. Let 
y = yv + yw, where W G E n ^ r (q m , n) such that V W = G¥(q m ) n . We have y v + (c 2 - y) v = c 2 , 
with rk(yy) < rk(y) < s and rk((c 2 — y)v) < rk(c 2 — y) < r — s. Therefore, rk(yy) = rk(y) = s, 
rk((c 2 - y) v ) = rk(c 2 - y) = r - s, and ©(yv) n ©((c 2 - y)y) = {0}. Since rk(y v ) = rk(y), we 
have ©(yw) Q ©(yv); and similarly ©((c 2 — y)w) Q ©(( c 2 — y)v)- Altogether, we obtain ©(yw) H 
©((c 2 - y)w) = {0}. However, yw + (c 2 - y)w = 0, and hence y w = (c 2 - y)w = 0. Therefore, 
y G V. 

We now prove that y is necessarily the projection of c 2 onto some ELS A of V. If y G V satisfies 
rk(y) = s and rk(c 2 — y ) = r — s, then y belongs to some ELS A and c 2 — y G B such that A © B = V. 
We hence have y = c 2 _4 and c 2 — y = c 2; #- 

On the other hand, for any A G E s (q m ,n) and £> G E r ^ s (q m ,n) such that .4 © i3 = V, c 2j ^ is a 
vector of rank s with distance r — s from c 2 by Lemma [3j By Lemma |U all the c 2 ^4 vectors are distinct. 
There are thus as many vectors y as ordered pairs (A,B). By Lemma [2l there are q s ( r ~ s ^ PI such pairs, 
and hence g s ( r ~ s ) [ r ] vectors y. ■ 

The problem of the intersection of three balls with rank radii is more complicated since the volume 
of the intersection of three balls with rank radii is not completely determined by the pairwise distances 
between the centers. We give a simple example to illustrate this point: consider GF(2 2 ) 3 and the vectors 



May 8, 2008 



DRAFT 



9 



ci = c'i = (0,0,0), C2 = c' 2 = (l,a,0), C3 = (a, 0,1), and c 3 = (a, a + 1,0), where a is a primitive 
element of the field. It can be verified that d R (ci,C2) = d R (c2,C3) = d R (c3,ci) = 2 and d R (c^,c 2 ) = 
4(c 2 ,c 3 ) = dj.(c' 3 ,c[) = 2. However, ^(ci) n £i(c 2 ) n Bi(c 3 ) = {(a + 1,0,0)}, whereas -B^c'J n 
-Bi(c' 2 ) n Bi(c' 3 ) = {(1, 0, 0), (0, a + 1, 0), (a, a, 0)}. We remark that this is similar to the problem of 
the intersection of three balls with Hamming radii discussed in [23, p. 58], provided that the underlying 
field is not GF(2). 

V. Covering properties of rank metric codes 
A. The sphere covering problem 

We denote the minimum cardinality of a code of length n and rank covering radius p as K R (q m , n,p). 
We remark that if C is a code over GF(q m ) with length n and covering radius p, then its transpose 
code C T is a code over GF(q n ) with length m and the same covering radius. Therefore, K R (q m , n, p) = 
K R (q n ,m, p), and without loss of generality we shall assume n < m henceforth in this section. Also 
note that K R (q m ,n,0) = q mn and K R (q m ,n,n) = 1 for all m and n. Hence we assume < p < n 
throughout this section. Two bounds on K R (q m ,n, p) can be easily derived. 

Proposition 5: For < p < n < m, < K R (q m ,n, p) < q m ( n ~p\ 

Proof: The lower bound is a straightforward generalization of the bound given in [10]. Note that 
the only codes with cardinality 2_ are perfect codes. However, there are no nontrivial perfect codes for 
the rank metric [8]. Therefore, K R (q m ,n, p) > Vjy. The upper bound follows from p < n — k for any 
(n, k) linear code (see [23] for a proof in the Hamming metric), and hence any linear code with covering 
radius p has cardinality < q m \ n ~P>. ■ 

We refer to the lower bound in Proposition [5] as the sphere covering bound. 

For a code over GF(q m ) with length n and covering radius < p < n, we have K R (q m ,n, p) < 
K H (q m , n, p), where K H (q m ,n, p) is the minimum cardinality of a (linear or nonlinear) code over GF(g m ) 
with length n and Hamming covering radius p. This holds because any code with Hamming covering 
radius p has rank covering radius < p. Since K H (q m ,n, p) < q m ( n ~P) [23], this provides a tighter bound 
than the one given in Proposition [5] 

Proposition 6: For < p < n < m, K R (q m ,n, p) > 3. 

Proof: Suppose there exists a code C of cardinality 2 and length n over GF(q m ) with covering 
radius p < n. Without loss of generality, we assume C = {0,c}. Since \B p (0) U B p (c)\ is a non- 
decreasing function of rk(c) by Corollary [T] we assume rk(c) = n. The code Q = (c) is hence an 
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(n, 1, n) lineal - MRD code over GF(g m ). Therefore, any codeword in Q\C is at distance n from C. Thus 
p = n, which contradicts our assumption. ■ 

Lemma 7: K R (q m ,n + n', p + p') < K R (q m , n, p)K R (q m , n', p') for all m > and nonnegative n, n' , 
p, and p'. In particular, K R (q m ,n + 1, p + 1) < K R (q m , n, p) and K R (q m , n + l,p) < q m K R (q m , n, p). 

Proof: For all x,y G G¥{q m ) n and x',y' G GF(g m ) n ', we have d R ((x,x'), (y,y')) < 4(x,y) + 
d R (y,y'). Therefore, for any C G GF(g m ) n , C" G GF(g m ) n ', we have p(C C") < p{C) + /9(C) and 
the first claim follows. In particular, (n',p') = (1,1) and (n',p') = (1,0) yield the other two claims 
respectively. ■ 

B. Lower bounds for the sphere covering problem 

We now derive two nontrivial lower bounds on K R (q m ,n, p). For < d < n, we denote the volume 
of the intersection of two balls in GF(q m ) n with rank radii p and a distance d between their respective 
centers d as I(q m ,n, p,d). When there is no ambiguity about the vector space considered, we simply 
denote it as I(p,d). I{p,d) is well defined by Lemma [6j and obviously I(p,d) = when d > 2p. 

Proposition 7: For < p < n < m and < I < \\og qm K R (q m ,n, p)\, 

M ,n ' P) " v{p)-I{p,n-l) • (1) 

Proof: Let us denote [logg™ K R {q m , n, p)\ as A for convenience. Let C = {cjj^Q 1 be a code of 

length n and covering radius p over GF(q m ). Define Cj = {cj^ =0 for < j < K - 1. For 1 < a < A 

and q m ( a -V < j < g ma , we have d R (cj,Cj-i) < n — a + 1 by the Singleton bound. The codeword 

Cj hence covers at most v{p) — I(p,n — a + 1) vectors that are not previously covered by Cj-\. For 

1 < / < A, the number of vectors covered by C thus satisfies 

q mn < v{p) + ^( ? ™ _ q(-^)[ v (p) - I(p, n-a + l)] + (K- q lm )Hp) ~ Hp, n-l)]. (2) 

a=l 

Since I(p, n — a + 1) = for a < n — 2p, © reduces to (Q]). ■ 
Note that the RHS of <[T]) is a non-decreasing function of /, thus the bound is tightest when I = 
[log qm K R (q m ,n, p)\. We obtain a lower bound by using the largest / such that the RHS of (fl} is less 
than 

Corollary 2: For < p < n < m, K R (q m ,n, p) > ^"^j . 

Proof: This is a special case of Proposition [7] for I = 0. ■ 
Corollary 3: For all m, n, and < p < [n/2\, K R (q m ,n, p) > — ^ ^ p . 

Proof: Since the balls of rank radius p around the codewords of a code with minimum rank 
distance 2p + 1 do not intersect, we have K R (q m ,n, p) > A R (q m , n, 2p + 1) = q m ( n ~ 2 P), and hence 
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2p > 0, and I(p, 2p) = q? 2 [ 



,r ' 2 p] by 



log„ m K R (q m , re, p) > n — 2p. Use Proposition [7] for I = n 
Proposition [4] ■ 
The bound in Corollary [3] can be viewed as the counterpart in the rank metric of the bound in [27, 
Theorem 1]. 

Van Wee [28], [29] derived several bounds on codes with Hamming covering radii based on the excess 
of a code, which is determined by the number of codewords covering the same vectors. Although the 
concepts in [28], [29] were developed for the Hamming metric, they are in fact independent of the 
underlying metric and thus are applicable to the rank metric as well. For all V C GF(q m ) n and a code 
C with covering radius p, the excess on V by C is defined to be Ec{V) = J2 c eC \^p( c ) H F| — |V|. 
If {Wi} is a family of disjoint subsets of GF(q m ) n , then E c (Ui^i) = T,i E c(Wi). Suppose Z = 
{z G GF(q m ) n \Ec{{z}) > 0}, i.e., Z is the set of vectors covered by at least two codewords in C. 
Note that z G Z if and only if \B p (z) n C\ > 2. It can be shown that \Z\ < E C {Z) = E c {GF{q m ) n ) = 
\C\V p {q m ,n) -q mn . 

Before deriving the second nontrivial lower bound, we need the following adaptation of [29, Lemma 

def 

8]. Let C be a code with length n and rank covering radius p over GF(q m ). We define A = {x G 
GF(<T)"|i(x } C) = p}. 

Lemma 8: For x G yl\Z and < p < n, we have that Ec{B\(~x)) > e, where 



def 



9" ft 1 " 1 



1 



+ (<f 



Proof: Since x ^ there is a unique Co G C such that d R (x, Co) = p. By Proposition [3] we have 
[S p (co)nBi(x)| = l + (g m -g p )[^]+(g p -l)[™]. For any codeword ci G C satisfying d R (x, ci) = p+1, 
by Proposition |4] we have |5 p (ci ) n B\ (x) [ = q p p^ 1 ] . Finally, for all other codewords c 2 G C at distance 

def 

> p + 1 from x, we have \B p (c2) n -Bi(x)| = 0. Denoting JV = |{ci G C|d R (x, c\) = p + 1}|, we obtain 
S c (Si(x)) = ^|B p (c)n J B 1 (x)|-| J B 1 (x)| 



cGC 



P+1 
1 



(<T " q") 



mod 



P+1 
1 



The proof is completed by realizing that (q m — q p ) ([?] — ["])< 0, while Ec{Bi{x)) is a non-negative 
integer. ■ 
For p = n — 1, Lemma [8] is improved to: 

Corollary 4: For x G A\2T and p = n-1, # c (.Bi(x)) = 0, where = f q n ~ l [™] (Cl-g™" 1 (<f" + [ n ^]). 
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Proof: The proof calls the same arguments as the proof above, with N = \C\ — 1 for p = n — 1. ■ 

Proposition 8: If e > 0, then K R (q m , n, p) > ^pjj^ n y where 5 = v(l) - qP' 1 [$] - 1 + 2e - 

The proof of Proposition [U provided in Appendix [A] uses the approach in the proof of [29, Theorem 
6] and is based on the concept of excess reviewed in Section III-BI The lower bounds in Propositions [7] 
and[U when applicable, are at least as tight as the sphere covering bound. For p = n — 1, Proposition [8] 
is refined into the following. 

Corollary 5: Let us denote the coefficients q n ~ 1 ["] and q n ~ l (g m + [ n ^ X ] ) as a and /3, respectively. 
K R (q m , n, re— 1) satisfies aK R (q m ,n, n—l) 2 —bK R (q m , n, n— l)+c > 0, where a = f a[v(n— l)+v(n— 2)], 
b = v(n-l) {q n ~ 2 - v(l) +p + l}+2aq mn +/3v(n-2), and c = q mn {2(3+l+q n - 2 f^ 1 ] -v(l)}. 

The proof of Corollary [5] is given in Appendix |B] 

C. Upper bounds for the sphere covering problem 

From the perspective of covering, the following lemma gives a characterization of MRD codes in terms 
of ELS's. 

Lemma 9: Let C be an (n, k) linear code over GF(q m ) (n < m). C is an MRD code if and only if 
C V = GF(q m ) n for all V G E n _ k {q m , n). 

Proof: Suppose C is an (n, k,n — k + 1) linear MRD code. It is clear that C n V = {0} and 
hence C V = GF(g m ) n for all V G E n _ k (q m ,n). Conversely, suppose C © V = GF(g m )" for all 
V G E n _i,(q m ,n). Then C does not contain any nonzero codeword of weight < n — k, and hence its 
minimum distance is n — k + 1. ■ 

For 1 < u < p, let a = 1, a x , . . . , a m+u -i G GF(q m+u ) be a basis set of GF(q m+u ) over GF(q), 
and let /3o = 1, . . . , /3 m _i be a basis of GF((7 m ) over GF(q). We define the linear mapping / between 

def 

two vector spaces GF(q m ) and & m = &(ao, ot\, . . . , a m _i) given by /(A) = «j for < i < m — 1. We 
remark that ao = 0o = 1 implies that / maps GF(g) to itself. This can be generalized to n-dimensional 
vectors, by applying / componentwise. We thus define / : GF(q m ) n — ► GF(q m+u ) n such that for 
any v = (vq, . . . , i> n -i), /(v) = (f(vo),...,f(v n -i)). Note that / depends on n, but we omit this 
dependence for simplicity of notation. This function / is a linear bijection from GF(q m ) n to its image 
S^, and hence / preserves the rank. / also introduces a connection between ELS's as shown below. 

Lemma 10: For u > 1, r < n, and any V G E r (q m ,n), /(V) C W, where W G E r (q rn+U ,n). 
Furthermore, / induces a bijection between E r (q m ,n) and E r (q m+U ,n). 

Proof: Let B = {bj} be an elementary basis of V G E r (q m , n). Then, bj G GF(g) n and bj = /(bj). 
Thus, {/(bj)} form an elementary basis, and hence /(V) C W, where W G E r (q m+U ,n) with {/(b,)} 



May 8, 2008 



DRAFT 



13 



as a basis. It is easy to verify that / induces a bijection between E r (q m ,n) and E r (q m+U ,n). ■ 
Proposition 9: Let C be an (n, n — p, p + 1) MRD code over GF{q m ) (n < m) with covering radius 
p. For < u < p, the code /(C), where / is as defined above, is a code of length n over GF(q m+u ) 
with cardinality q m ( n ~P^ and covering radius p. 

Proof: The other parameters for the code are obvious, and it suffices to establish the covering 
radius. Let X u be a subspace of GF(q m+u ) with dimension u such that 6 m 1 U = GF(q m+u ). Any 
u G GF(q m+u ) n can be expressed as u = v + w, where v G SJ^ and w G X™. Hence rk(w) < u, 
and w G W for some W G E p (q m+U ,n) by Lemma Q] By Lemmas [9] and [lOl we can express v as 
v = /(c + e) = /(c) + /(e), where c G C and e G V, such that /(V) C W. Eventually, we have 
u = /(c) + /(e) + w, where /(e) + w G W, and thus d(u, /(c)) < p. Thus /(C) has covering radius 
< p. Finally, it is easy to verify that the covering radius of / (C) is exactly p. ■ 
Corollary 6: For < p < n < m, K R (q m ,n,p) < g max{m- P ,n}(n- P )_ 

Proof: We can construct an (n,n — p) linear MRD code C over GF(q^) with covering radius 
p, where p = max{m — p, n}. By Proposition |9j f(C) C GF(g m ) n , where / is a rank-preserving 
mapping from GF(q fl ) n to a subset of GF(q m ) n similar to / above, has covering radius < p. Thus, 
K R (q m ,n,p) < |/(0| = \C\=q^ n -P\ ■ 
We use the properties of K R (q m ,n, p) in Lemma [7] in order to obtain a tighter upper bound when 
p > m — n. 

Proposition 10: Given fixed m, n, and p, for any < I < n and (rii,pi) for < i < I — 1 so that 
< rii < n, < pi < rii, and nj + pi <m for all i, and X^=o n-i = n and X^to Z 5 * = Z 9 ' we nave 

^T R (g m ,n,p)< min ( 3 ) 
Proof: By Lemma |7J we have K R (q m ,n, p) < Y\ i K R (q m ,ni,pi) for all possible sequences {^j 
and {nj}. For all i, we have K R (q m ,rii, pi) < q^-P^-P*) by Corollary and hence K R (q m ,n, p) < 

It is clear that the upper bound in ® is tighter than the upper bound in Proposition [5] It can also be 
shown that it is tighter than the bound in Corollary [6] 

The following upper bound is an adaptation of [23, Theorem 12.1.2]. 
Proposition 11: For < p < n < m, K R (q m , n, p) < !_ log (gmn _^ p)) + L 

Our proof, given in Appendix ICl adopts the approach used to prove [23, Theorem 12.1.2]. Refining 
[23, Theorem 12.2.1] for the rank metric, we obtain the following upper bound. 
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def 

Proposition 12: For < p < n < m and a = min{n, 2p}, 



,11171 



K R (q m ,n,p) < k v() ——^ — -- — + ——H min{s j} , (4) 

w \min{s,jj v{pj J v(p) 1 ,J/ 

where k v(p) = q mn - v(p)q m{ - n - a \ j = \v(p) - v(p) 2 q~ ma ^, s = v(p) - £f =(l _ p g i(a_i) [;]> and H k = 
2~Zi=i j is the k-th harmonic number. 

Proof: We denote the vectors of GF(q m ) n as Vj for i = 0, 1, ... , q mn — 1 and we consider a 
q mn x q mn square matrix A defined as ay = 1 if d R (vj, Vj) < p and ay = otherwise. Note that each 
row and each column of A has exactly v(p) ones. We present an algorithm that selects K columns of A 
with no all-zero rows. These columns thus represent a code with cardinality K and covering radius p. 

Set A„( p ) = A and k v ^ + i = q mn . For i = v{p),v{p) — 1,... ,1, repeat the following step: First, 
select from Aj a maximal set of Ki columns of weight i with pairwise disjoint supports; Then, remove 
these columns and all the ii-Q = ki + \ — ki rows incident to one of them, and denote the remaining 
ki x (qr mn — K v i p \ — ■ ■ ■ — Ki) matrix as Aj_i. The set of all selected columns hence contains no all-zero 
rows. 

For 2p < n, we can select an (n, n — 2p, 2p + 1) linear MRD code C for the first step. Since MRD 
codes are maximal, they satisfy the condition. If 2p > n, C is chosen to be a single codeword, and 
K v(p) = 1. Thus K v(p) = q m ( n - a ) and k v{p) = q mn - v(p)q m ( n ~ a \ where a = min{rt,2p}. 

We now establish two upper bounds on ki for 1 < i < v(p). First, it is obvious that ki < k v (p). Also, 
every row of Aj_i contains exactly v(p) ones; on the other hand, every column of Aj_i contains at most 
% — 1 ones. Hence for 1 < i < v{p), v{p)ki < (i — l)(q mn — K v ^ — ■ ■ ■ — Ki) < (i — l)q mn , and thus 



..inn 



h<(i-i)^. (5) 

v(p) 



def 



v(p)k v{pj 



Clearly k x = by ©. We have k v{p) < (i - 1)^ if i - 1 > j 

We now establish an upper bound on K = Y!,i=i Ki. For any vector x £ GF(g m ) n , we have d R (x, C) < 
2p. This is trivial for 2p > n, while for 2p < n this is because MRD codes are maximal codes. For 2p < n, 
at least q p2 [ 2 ^] vectors in B p (x) are already covered by C by Proposition [4] For 2p > n, it can be shown 
that at least 2~2i= n -p ( l^ n [1] vectors in B p (x) are already covered by C. It follows that after the first 
step, the column weight is at most s = f v(p)— Y^i= a -p q^ a ~ 1 ^ ["] ■ Since s > max{l < i < v{p) : K t > 0}, 
K = K v r p) + Yfi=\ K i = K v(p) + ^r 1 + 2~2t=2 Kftj\- Usin § the two upper bounds on ki above, we 

Obtain K < k v{p) - ife) + ^ H min{s,j}. ■ 

Following [23], where [23, Theorem 12.1.2] is referred to as the Johnson-Stein-Lovasz theorem, 
we refer to the algorithm described in the proof of Proposition [12] as the Johnson-Stein-Lovasz (JSL) 
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algorithm. The upper bound in Proposition [12] can be loosened into the following. 

Corollary 7: For < p < n < m, K R (q m ,n, p) < gy (lnv(p) + 7+ 2 v(p)+i/3 ) ' wnere 7 * s tne 
Euler-Mascheroni constant [30]. 

Proof: Using ©, we obtain K = ^ < gl + ijfey < ^H v{p) . The proof is 

concluded by H v{p) < lnv(p) + 7 + 2v {p)+i/z t 3 °]- ■ 

D. Covering properties of linear rank metric codes 

Proposition [5] yields bounds on the dimension of a linear code with a given rank covering radius. 
Proposition 13: An (n,k) linear code over GF(q m ) with rank covering radius p satisfies n — p — 

p( n- P )-lo gq K 

m — " 

Proof: The upper bound directly follows the upper bound in Proposition [5] We now prove the 
lower bound. By the sphere covering bound, we have q mk > However, by Lemma [5] we have 
v (p) < q p( - m+n ~^~ log i Kq and hence q mk > q mn -p( m + n -p)+ lo s q K i_ B 

We do not adapt the bounds in Propositions [7] and [8] as their advantage over the lower bound in 
Proposition [13] is not significant. Next, we show that the dimension of a linear code with a given rank 
covering radius can be determined under some conditions. 

Proposition 14: Let C be an (n, k) linear code over GF(g m ) (n < m) with rank covering radius p. 
Then k = n — p if p G {0, 1, n — 1, n} or p(n — p) < m + \og q K q , or if C is a generalized Gabidulin 
code or an ELS. 

Proof: The cases p G {0, n — l,n} are straightforward. In all other cases, since k < n — p by 
Proposition [13] it suffices to prove that k > n — p. First, suppose p = 1, then fc satisfies <7 mfc > 
by the sphere covering bound. However, v(l) < q m+n < q 2m , and hence k > n — 2. Second, if 
p(n — p) < m + logg i^g, then < ^ (p(n — p) — \og q K q ^ < 1 and k > n — p by Proposition [T3l Third, 
if C is an (n,k,n — k + 1) generalized Gabidulin code with k < n, then there exists an (n, k + l,n — k) 
generalized Gabidulin code C such that C C C. We have p > d R (C) = n — k, as noted in Section ITl-B I 
and hence k > n — p. The case k = n is straightforward. Finally, if C is an ELS of dimension k, then 
for all x with rank n and for any c G C, d R (x, c) > rk(x) — rk(c) > n — k. ■ 
A similar argument can be used to bound the covering radius of the cartesian products of generalized 
Gabidulin codes. 

Corollary 8: Let Q be an (n,k,d R ) generalized Gabidulin code (n < m), and let Q l be the code 
obtained by I cartesian products of Q for / > 1. Then the rank covering radius of Q l satisfies p{G l ) > d R —l. 
Note that when n = m, Q l is a maximal code, and hence Corollary [8] can be further strengthened. 
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Corollary 9: Let Q be an (m, fc, g? r ) generalized Gabidulin code over GF(g m ), and let Q l be the code 
obtained by I cartesian products of Q. Then p(Q ) = d R — 1. 

E Numerical methods 

In addition to the above bounds, we use several different numerical methods to obtain tighter upper 
bounds for relatively small values of m, n, and p. First, the JSL algorithm described in the proof of 
Proposition [12] is implemented for small parameter values. Second, local search algorithms [23] similar 
to the ones available for Hamming metric codes are somewhat less complex than the JSL algorithm. 
Although the complexity for large parameter values is prohibitive, it is feasible. Third, we construct 
linear codes with good covering properties, because linear codes have lower complexity. 

We can verify if a covering radius is achievable by a given code size by brute force verification, thereby 
establishing lower bounds on K R (q m ,n, p). Obviously, this is practical for only small parameter values. 

F. Tables 

In TableJl we provide bounds on K R (q m , n, p), for 2 < m < 7, 2 < n < m, and 1 < p < 6. Obviously, 
K R (q m ,n, p) = 1 when p = n. For other sets of parameters, the tightest lower and upper bounds on 
K R (q m ,n, p) are given, and letters associated with the numbers are used to indicate the tightest bound. 
The lower case letters a-f correspond to the lower bounds in Propositions [5J EJ and [71 Corollaries [2] and [3j 
and Proposition [8] respectively. The lower case letter g corresponds to lower bounds obtained by brute 
force verification. The upper case letters A-E denote the upper bounds in Proposition [5J Corollary [6] 
and Propositions [10J [TT] and [12] respectively. The upper case letters F-H correspond to upper bounds 
obtained by the JSL algorithm, local search algorithm, and explicit linear constructions respectively. 

In Table [III we provide bounds on the minimum dimension k for q = 2, 4 < m < 8, 4 < n < m, and 
2 < p < 6. The unmarked entries correspond to Proposition [14] The lower case letters a and e correspond 
to the lower bound in Proposition [13] and the adaptation of Corollary [3j to linear codes respectively. The 
lower case letter h corresponds to lower bounds obtained by brute force verification for linear codes. The 
upper case letter A corresponds to the upper bound in Proposition [13] The upper case letter H corresponds 
to upper bounds obtained by explicit linear constructions. 

Although no analytical expression for I(p, d) is known to us, it can be obtained by simple counting 
for the bounds in Proposition [7J or Corollary [3] In Appendix |Pl we present the values of I(q m ,n, p,d) 
used in calculating the values of the bounds in Proposition [7J and Corollary [2] displayed in Table [U We 



May 8, 2008 



DRAFT 



17 



171 


71 


„ 1 

P — 1 


p — z 


p — 6 


„ A 

p — 4 


p — b 


P 


— 


2 


2 


e o r 


i 
l 












3 


2 


e 4 d 


t 
1 














J 


a i 1 i £ n 
e 1 l-lo r 


g 4 C 


i 
1 










A 

4 


2 


.-.TOD 

e /-o rS 


t 
1 














3 


e 4U-64 b 


d 4-7 F 


l 












4 


f om Tin t? 
I Ay 5- 1 12. r 


e 10-48 (j 


D 5-1 r 


i 










2 


e Iz-lo n 


t 
1 














3 


e 154-z5o B 


d 6-8 B 


i 
1 












4 


e zzo/-4U9o rs 


e 33-256 (_ 


d 4-8 C 


i 
i 










^ 

j 


e 34oV4-z L. 


toq Tool c 

e 2jJ-2S81 b 


„ n qo 1 1 
a y- 51 ti 


b j-8 C 


i 






6 


2 


e 23-32 B 


1 














3 


e oUl-ll)z4 d 


d 1 1-16 B 


1 












4 


„ 1 "7QOO o!5 D 

e 1 /ozz-z d 


C 124-256 -b 


a o-lo L. 


t 
1 










5 


e 550395-2 B 


e 1770-2 C 


1 31-256 C 


a 3-16 C 


i 








6 


I 1 /318410-2 L 


I 2/065-413582 b 


1 214-421 1 b 


f 9-181 D 


b 3-16 C 




1 


7 


2 


e 44-64 B 


1 














3 


e 23 / 2-4096 B 


d 20-32 B 


1 












4 


e 141231-2 18 B 


f 484-1024 B 


a 9-16 B 


1 










5 


e 8735289-2 24 B 


e 13835-2 15 B 


a 111-1024 C 


a 5-16 C 


1 








6 


e 549829402-2 30 B 


f 42229-2 22 C 


e 1584-2 15 C 


f 29-734 E 


a 3-16 C 




1 




7 


e 34901004402-2 37 C 


f 13205450-239280759 E 


e 23978-586397 E 


f 203-5806 E 


a 8-242 D 


b 3 


■16 C 



TABLE I 

BOUNDS ON K R (q m , n, p), FOR 2 <m <7, 2 <n <m, AND 1 < p < 6. 



also present the codes, obtained by the numerical methods in Section IV-El that achieve the tightest upper 
bounds in Tables U and ITT] 

G. Asymptotic covering properties 

Table U provides solutions to the sphere covering problem for only small values of m, n, and p. Next, 
we study the asymptotic covering properties when both block length and minimum rank distance go to 
infinity. As in Section [TTlJ we consider the case where lim^^oo ^ = b, where b is a constant. In other 
words, these asymptotic covering properties provide insights on the covering properties of long rank 
metric codes over large fields. 
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~ /I C A 

a 4-j A 


a 


3-4 A 


a Z-j A 


~ 1 O A 

a 1-2 A 


1 


8 


4 


2 




1 











5 


3 




2 


1 









6 


a 3-4 A 


a 


2-3 A 


a 1-2 A 


1 







7 


a 4-5 A 


a 


3-4 A 


a 2-3 A 


a 1-2 A 


1 




8 


a 5-6 A 


a 


3-5 A 


a 2-4 A 


a 1-3 A 


a 1-2 A 



TABLE II 

Bounds on k for g = 2, 4 < m < 8, 4 < n < m, and 2 < p < 6. 



The asymptotic form of the bounds in Lemma [5] are given in the lemma below. 

Lemma 11: For < 5 < min{l, b^ 1 }, lim n ^oo [log 9m „ V\g n \ (q m , n)] = 5(1 + b — b5). 

Proof: By Lemma we have q Mm+n-d R ) < y ^ < if -i^ R (m+n-d R )_ Taking the logarithm, this 

becomes 5(1 + b — b5) < log gm „ v ([5n\) < 5(1 + b — b5) — — f^- 1 . The proof is concluded by taking 
the limit when n tends to infinity. ■ 

Define r = f ^ and A:(r) = linin^oo inf [log gm n K R (q m , n, p)] . The bounds in Proposition [5] and 
Corollary [7] together solve the asymptotic sphere covering problem. 

Theorem 1: For all b and r, A;(r) = (1 — r)(l — 6r). 

Proof: By Lemma [TTI the sphere covering bound asymptotically becomes k(r) > (1 — r)(l — br). 
Also, by Corollary|7J K R (q m ,n, p) < [1 + ln'u(p)] < [1 + mralng] and hence log qmn K R (q m , n, p) < 
log ?m „ + 0((mn) _1 lnmn). By Lemma [TTI this asymptotically becomes k(r) < (1 — r)(l — 6r). Note 
that although we assume n < m above for convenience, both bounds in Proposition [5] and Corollary |7] 
hold for any values of m and n. ■ 



May 8, 2008 



DRAFT 



19 



Appendix 

A. Proof of Proposition [S] 

We first establish a key lemma. 

Lemma 12: If z E Z and < p < n, then \A n #i(z)| < v(l) - qp- 1 [?] . 

Proof: By definition of p, there exists c G C such that d R (z,c) < p. By Proposition 12 |-£>i(z) n 
Bp_i(c)| gets its minimal value for d R (z, c) = p, which is q p ~ l [£] by Proposition [4] A vector at distance 
< p — 1 from any codeword does not belong to A. Therefore, B\{z) n B p -\(c) C Bi(z)\A, and hence 
|A n B 1 (z)\ = (z) | - |Bi(z)\A| < u(l) - |Bi(z) n Bp_i(c)|. ■ 
We now give a proof of Proposition [8] 

Proof: For a code C with covering radius p and e > 1, 

def 



7 = ^^.icr^.^.^.^^l^. 

< e|A|-(e-l)|Z| 

< e|A| - (e - l)\A<T\ Z\ = e\A\Z\ + \AnZ 
where $7} follows from \Z\ < \C\v(p) — q mn , given in Section Hl-Bl 



(6) 
(V) 



a£A\Z aeAnZ 
aGA 

where ([8]> follows from Lemma [8] and \Af]Z\ < E C {A n Z). 



(8) 



^ ^ E E ^(W) 

aGAxG-Bi(a)nZ 

= E E E c ({x}) = Y J \AnB 1 (x)\E c ({x}\ 

xGZaGBi(x)nA xGZ 

where (O follows the fact that the second summation is over disjoint sets {x}. By Lemma 



(9) 



we obtain 



7 < ErW-^ 1 



xez v 
v(l) - qP- 1 



v(l) — q 



p-i 



ljj£c({x}) 

E C {Z) 

(\C\v(p)-q« 



(10) 



Combining (fTOl and ©, we obtain the bound in Proposition [8] 
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B. Proof of Corollary \5\ 

For p = n - 1, CE!) becomes [q mn - \C\v(n - 2)] - ((f> - 1) [|C|u(n - 1) - <? mn ] < (u(l) - 
qn-2\n- Y)(\C\v(n — 1) — q mn ). Substituting (f> = a\C\ — j3 and rearranging, we obtain the quadratic 
inequality in Corollary [5] 

C. Proof of Proposition [77] 

Given a radius p and a code C, denote the set of vectors in GF(q m ) n at distance > p from C as 
P p (C). To simplify notations, Q = f q mn and p p (C) == <3 _1 |P p (C)|. Let us denote the set of all codes 
over GF(q m ) of length n and cardinality K as Sk- Clearly \Sk\ = ( K )- The average value of p p (C) 
for all codes C 6 S*x is given by 

- ^5J p ' (c)l ^nS„ e X,>/ 

1 A| x£F C&S K \d R (x,C)>p 

Q-*(p)\ //Q 

Eq. (fTTb comes from the fact that there are (®~k P ^) codes with cardinality if that do not cover x. For 
all K, there exists a code C £ Sk for which p p (C) is no more than the average, that is: 

Let us choose if = — log fTZ^pr^^ + 1 so that -^1°Sq (l ~~ Q~ lv (p)) < — 1 and hence p p {C) = 
(l — Q~ 1 v(p)) K < Q^ 1 . It follows that |P p (C")[ < 1, and C has covering radius at most p. 

D. Numerical results 

The values of I(q m ,n, p,d), used in calculating the bounds in Proposition UJ or Corollary [2j obtained 
by counting are i(2 4 ,3,2,3) = 560, i(2 5 ,3,2,3) = 1232, i(2 5 ,4,3,4) = 31040, i(2 6 ,3,2,3) = 2576, 
i(2 6 , 4, 2, 3) = 2912, i(2 6 , 4, 3, 4) = 756800, and i(2 7 , 3, 2, 3) = 5264. 

We now present the codes, obtained by computer search, that achieve the tightest upper bounds in 
Tables HI and ITT1 The finite fields use the default generator polynomials from MATLAB [31]. First, the 
linear code used to show that if R (2 5 ,5,3) < 32 has a generator matrix given by G = (1, a, a 2 , 0, 0), 
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where a is a primitive element of GF(2 5 ). We use the skip-vector form [32] to represent the other codes 
obtained by computer search. The skip-vector form of a code C = {cj}^ 1 over G¥(q m ) n can be 
obtained as follows. First, each codeword Cj € G¥{q m ) n is represented by an integer xi in [0, q mn — 1] 
according to the lexicographical order. Second, the integers x« are sorted in ascending order; the resulting 
integers are denoted as x\. Third, calculate yi defined as yo = x' and yi = x\ — x\_ x — 1 for 1 < % < K — l. 
Fourth, if yi = y i+1 = ... = y i+k -i, then we write y\. 

Below are the codes obtained by the JSL algorithm. 
K R (2 2 ,2,1)= 3 3 

K R (2 3 , 3, 1) < 16 5 25 66 21 51 9 20 5 85 21 2 25 49 9 84 5 
K R (2 4 , 3, 2) < 7 135 689 34 420 477 522 759 

K R (2 4 ,4, 1) < 722 57 308 349 86 125 27 192 18 38 36 95 86 64 157 67 98 7 301 21 131 42 
39 60 149 97 40 116 20 85 11 90 15 56 19 167 9 137 10 7 75 21 51 18 110 12 82 27 38 15 143 2 24 
120 5 39 77 223 14 52 27 12 179 42 86 88 100 3 130 34 66 35 5 30 171 210 137 34 29 149 59 69 97 
26 105 93 286 61 14 30 136 62 148 97 132 182 184 3 69 43 31 74 5 190 85 92 85 91 146 58 75 20 
4 234 292 56 10 40 56 86 37 85 111 80 32 103 69 82 34 106 187 42 44 47 242 220 43 10 144 27 50 
97 118 60 94 61 297 36 12 222 19 16 88 72 170 19 14 197 20 120 136 54 20 59 47 86 49 37 70 216 
164 2 92 53 77 83 70 225 73 38 119 33 224 34 316 1 51 36 74 33 19 128 60 52 160 31 62 135 50 135 
282 19 38 140 80 88 55 65 50 46 22 16 320 15 110 58 183 106 30 170 128 82 2 152 189 60 62 61 
180 30 74 22 15 201 16 184 44 206 59 93 16 148 12 94 33 102 40 68 52 12 114 32 216 45 134 31 
140 29 324 87 97 206 14 26 42 4 22 48 89 60 85 29 14 203 37 7 300 165 128 58 224 80 95 3 22 98 
90 4 337 6 25 121 64 54 84 13 109 87 30 49 32 56 26 116 40 126 109 47 27 100 68 14 98 60 167 33 
90 224 6 229 262 89 48 89 63 157 107 21 28 445 2 13 26 132 7 36 3 81 11 50 43 35 127 89 7 180 
26 22 89 82 18 113 230 49 278 197 323 24 93 230 144 99 15 8 255 27 9 19 79 80 56 175 107 40 62 
105 20 115 9 41 95 72 97 109 250 51 166 47 65 94 7 166 133 108 148 56 76 201 69 98 133 33 46 13 
36 176 12 44 20 23 90 96 98 191 56 90 162 66 39 44 107 198 90 124 353 354 242 21 170 161 35 
211 9 14 5 155 13 20 4 120 24 89 36 73 139 98 114 128 30 64 33 67 132 15 102 105 22 48 161 36 
35 53 19 102 150 4 30 54 18 119 14 19 60 84 2 50 62 40 95 13 33 140 38 28 116 60 167 44 104 
244 366 93 87 9 282 157 158 248 19 7 123 218 2 130 236 178 13 12 46 97 67 30 98 25 26 49 111 
40 36 197 2 58 67 18 98 155 21 34 9 93 101 61 8 111 71 68 112 232 69 403 9 148 40 237 248 99 
93 230 53 171 49 89 131 13 110 27 157 107 58 19 16 19 92 110 366 68 81 198 212 73 57 193 158 
33 123 129 52 85 23 181 48 85 150 200 73 74 41 36 183 79 72 278 145 240 26 27 144 49 212 99 82 
173 93 221 118 33 108 39 ll 2 122 20 4 12 136 177 45 39 51 6 150 30 50 10 228 4 146 77 14 78 
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117 88 141 39 260 358 1 97 170 39 248 116 30 118 15 11 49 271 83 8 118 32 54 96 21 67 71 234 97 
229 106 59 166 19 35 152 42 56 317 11 184 90 2 60 65 15 272 231 121 56 53 11 93 250 272 38 26 
88 6 110 59 158 14 109 29 110 113 58 206 87 46 162 99 13 22 59 220 146 161 152 73 69 162 

Below is the code obtained by a local search algorithm. 
K R (2 4 ,4,2) < 48 1493 1124 265 285 1030 2524 1366 493 6079 968 2145 848 312 473 1307 712 
1088 2274 1380 1114 1028 567 422 1462 699 203 180 4669 146 978 3933 1810 2083 345 354 659 
1054 2314 1443 2660 2675 1512 756 1229 95 2144 1624 1148 
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